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ABSTRACT 

We present results of the application of turbulence models and 
the Reynolds analogy to the Navier-Stokes computations of Mach 
2.9 two-dimensional compression ramp flows. We studied the 
Baldwin-Lomax eddy viscosity model and the k — E turbulence 
transport equations for the turbulent momentum flux modeling in 
the Navier-Stokes equations. We also studied the Reynolds analogy 
for the turbulent heat flux modeling in the energy equation. The 
Navier-Stokes equations and the energy equation were numerically 
solved for the flow properties. The Reynolds shear stress, the skin 
friction factor, and the surface heat transfer rate were calculated 
and compared with their measurements. We concluded that (a) with 
a hybrid k—E turbulence model for turbulence modeling, the 
present computations predicted the skin friction factors of the 8° 
and 16° compression ramp flows and (b) with the turbulent Prandtl 
number Pr, = 0.93 and the ratio of the turbulent thermal and mo- 
mentum transport coefficients \^q^t = 2/Pr^ the present com- 
putations also predicted the surface heat transfer rates beneath the 
boundary layer flow of the 16° compression ramp. 

NOMENCLATURES 

Cy skin friction factor, T w / (1 / 2)p x U] o 

C i’ C tf’ C klcb’ C g’ C j sempirical constants 

C p specific heat 

K thermal conductivity 

k turbulence kinetic energy 

n direction normal to wall surface 


p 

static pressure 


laminar Prandtl number, 0.73 


turbulent Prandtl number, 0.93 

Q 

thermal flux 

Re 

Reynolds number, p x U„5 0 f p„ 

S 

Reynolds analogy factor, {^p q l Pr, 

T 

static temperature 

t 

time 

t r 

temperature fluctuation 

U 

X-direction velocity component 


/ \ 0 - 5 
frictional velocity, \Ji w \(0\ w f p w j 

u. 

velocity component parallel to wall 
surface 

u f 

X-direction velocity fluctuation 

V 

T-direction velocity component 

v' 

T-direction velocity fluctuation 

X,Y 

physical coordinates 
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boundary layer thickness 
Kronecker delta 
turbulence dissipation rate 
dynamic viscosity 

turbulent thermal transport coefficient 

turbulent eddy viscosity 

density 

surface shear stress 

Reynolds shear stress, f^ [ (dU / dY + <?V / dX) 
vorticity 

condition at ramp comer 
downstream of ramp comer 
sublayer condition 
upstream of ramp comer 
sublayer edge condition 
wall surface condition 
condition at X = 0 


and numerical solutions of the compressible Navier-Stokes equa- 
tions. Existing computational studies [2 and 3] indicate that, with 
proper turbulence models, solutions of the Navier-Stokes equations 
and the energy equation could predict the surface pressure and skin 
friction variations for most of the interaction flows. 

The surface heat transfer within the shock wave/boundary layer 
interaction flow is also very important to the design of a high speed 
propulsion inlet system. For many turbulent flow analyses, a prac- 
tical approach to estimate the turbulence effect on the surface heat 
transfer has been the Reynolds analogy. However, boundary layer 
theory [4] indicates that the use of Reynolds analogy requires a 
wide range of turbulent Prandtl numbers in order to properly model 
the turbulent thermal fluxes in pipe flow and a heated boundary 
layer flow. Closer to the present work, a recent experiment [5] of a 
supersonic compression ramp flow also revealed that the Reynolds 
analogy needed an additional factor to correlate the Stanton num- 
ber and the local skin friction factor Therefore, there is a need to 
investigate the application of the Reynolds analogy for the model- 
ing of the turbulent thermal flux within a shock wave/turbulent 
boundary layer interaction flow field. 

The present paper describes some Navier-Stokes computational 
analyses of the Reynolds shear stresses and the surface heat trans- 
fer in a two-dimensional shock wave/turbulent boundary layer in- 
teraction flow. We considered the shock wave/boundary layer inter- 
action around 8° and 16° compression ramps (Fig. 1 ) imbedded in a 
fully developed turbulent boundary layer flow. The free stream Mach 
number was 2.9. The free stream Reynolds number (based on the 
initial boundary layer thickness & 0 ) was 1.7xl0 6 . These are the ex- 
perimental conditions of Ref. [5]. The Baldwin-Lomax turbulence 
model [6], the low Reynolds number k - £ turbulence model of 
Nichols [71 and their revised versions [8] were each used, respec- 
tively, to model the turbulent momentum fluxes in the momentum 
equation. The Reynolds analogy was then used to model the turbu- 
lent heat flux in the energy equation. The turbulent Prandtl number 
was constant and equal to 0.93. The turbulent thermal flux was rep- 
resented by the product of the mean temperature gradient and a 
transport coefficient. A Reynolds analogy factor was used to ac- 
count for the stream-wise variation of the ratio of the thermal trans- 


oo free stream condition 

0,00 free stream stagnation condition 


incompressible flow condition 

INTRODUCTION 

Shock wave/turbulent boundary layer interaction flow is an im- 
portant phenomenon which occurs in the flow field of a high speed 
propulsion inlet system. A review of the existing literature on the 
shock wave/boundary layer interaction flows can be found in an 
AGARD report [1]. Generally, the existing studies have described 
the flow phenomena for transonic and supersonic speed regions by 
using experimental observations, correlations, theoretical concepts, 
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Figure 1 . — Compression ramp flow model (8 0 = 26 mm, 
Re = Poo U«VMoo = 1*76x106). 
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port coefficient and the turbulent eddy viscosity. The compressible 
Navier-Stokes equations and the energy equation were solved for 
the flow properties with a time-dependent explicit finite difference 
computational code [9], 

The present Navier-Stokes computational results of the turbulent 
eddy viscosity, the Reynolds shear stress, the mean flow velocity, 
and the surface heat transfer rate within the interaction flow fields 
are described and compared with their measurements in the exist- 
ing experiments [5, 10, and 11]. Based on the comparisons, the ef- 
fectiveness of the turbulence modeling techniques and the Reynolds 
analogy on the predictions of the Reynolds shear stress, the skin 
friction factor, the surface pressure, and the surface heat transfer 
rate of the supersonic compression ramp flow are discussed in the 
paper. 

TURBULENCE MODELS 

The Baldwin and Lomax model [6] is widely used for turbulence 
modeling in compressible flows. A compressible k — £ two equa- 
tion turbulence model developed by Nichols has been used in some 
supersonic flow computations [7]. Viscous sublayer turbulence 
models [12 and 13] have been developed and incorporated with the 
k and e equations to resolve the near wall turbulence modeling. 
Previous studies [7, 12, and 13] presented the applications of the 
k - e turbulence model and the viscous sublayer turbulence mod- 
els in the predictions of the surface pressure and skin frictions of 
compression ramp flows with large deflection angle and free stream 
Mach number. The effectiveness of these models on the computa- 
tions of the Reynolds shear stress and the surface heat transfer rate 
within a compression ramp flow were not reported. 

In the present work, we studied the existing turbulence modeling 
techniques [6, 7, 12, and 13] for the computations of the Reynolds 
shear stresses and surface heat transfer rates of two-dimensional 
compression ramp flows (Fig. 1). Two flow Fields with different 
deflection angles (8° and 16°) were considered. Experiments [10] 
of the flow fields indicated that flow separation did not occur. Thus, 
we did not consider the turbulence modeling of a separated flow. A 
summary of the turbulence modeling techniques used here is given 
in the following section. 

Baldwin and Lomax Turbulence Model 

This turbulence model is a two-layer eddy viscosity model. The 
eddy viscosity formulas of a two-dimensional flat plate turbulent 
boundary layer flow are shown in the following: 

In the inner layer, the eddy viscosity, |i r •, is described by 



In the outer layer, the eddy viscosity, p, is described by 


/!,„ =0.0168 C.pF^F^, (2) 


F wakc is thc smalIer of y max F m« and 025 Y m m U Lj F m m ' 


««.=(Vi'’+V ! )_ B ,and F=yH(l-e- r “) 

The quantity, F max , is the maximum value of F that occurs in a 
profile and F max is the value of Y at which F max occurs. The func- 
tion, F kleb , is the Klebanoff intermittency factor given by 

F “ = i +5 .5(c„r/r M ,)* < 3 > 

In earlier work [8], the empirical constants, C-, C Q , and C k)eb , 
appearing in Eqs. (1) to (3) were adjusted to achieve better predic- 
tion of the Reynolds shear stress in the supersonic compression ramp 
flow field. Different sets of these empirical constants were found to 
model the turbulence eddy viscosity upstream and downstream of 
the compression ramp. These sets of empirical constants are (1) 
Cy = 0.41, C Q - 1, and C Ueb = 0.8 for the upstream modeling, 
X < X c and (2) C i = 0.45, C G = 2, and = 0.8 for the downstream 
modeling, X>X C . Using these empirical constants, we modified the 
Clauser constant and the Klebanoff intermittency factor of the 
Baldwin-Lomax turbulence model. 

k - e Turbulence Models 

The following compressible low Reynolds number turbulence 
model was first used [7] to predict the turbulent eddy viscosity within 
the supersonic flat plate turbulent boundary layer flow. 

The turbulent eddy viscosity, p r , was defined as 

,4, 


with = 0.09. 

The turbulence kinetic energy, k , and the dissipation rate, £, were 
described by using the following transport equations: 
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The production term, P, was defined as 
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k v and e v are the turbulence kinetic energy and its dissipation rate 
at the viscous sublayer edge. They are 
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Equations (4) to (8) were also solved for the turbulent eddy vis- 
cosity in the computations of the interaction flow along the com- 
pression ramp. The source terms were found to be numerically stiff 
at locations near the comer of the compression ramp. Therefore, a 
viscous sublayer model [12] was used to eliminate the stiffness prob- 
lem. This viscous sublayer model and the high Reynolds number 
forms of the k — t turbulence model corresponding to Eqs. (4) to 
(8) were combined to formulate the following hybrid k~Z turbu- 
lence model [8] for the flow over the compression ramp. 

Adjacent to the wall, the viscous sublayer thickness, F v , was de- 
fined by 


Y v 



f £ 0.23 


u. 


(9) 


where C s (= 0.4) is an empirical constant. Within the viscous sublayer 
thickness, 0 < Y < Y v , the turbulent eddy viscosity, \i t s , was ex- 
pressed as 


At locations outside of the viscous sublayer, the high-Reynolds- 
number forms of Eqs. (4) to (8) were solved for the eddy viscosity. 
The free stream turbulence kinetic energy and the dissipation rate 
were the far field boundary conditions. Similar turbulence quanti- 
ties given by Eqs. (13) and (14) were the boundary conditions at the 
viscous sublayer edge. 

REYNOLDS ANALOGY 

The general and practical heat transfer modeling technique for a 
turbulent flow is the Reynolds analogy. Given a turbulent Prandtl 
number, the analogy relates the turbulent thermal flux to the turbu- 
lent momentum flux. In this section, we describe the application of 
the Reynolds analogy in the prediction of the surface heat transfer 

downstream of the compression ramp flow. _ 

We postulate that the turbulent thermal fluxes, pu't' and pv't' 
assume the forms 
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where the turbulence kinetic energy, k s , and the dissipation rate, 
e,.are 



( 11 ) 


pvr=_p *ar (16) 

where \i q is the thermal transport coefficient. Including the thermal 
conduction fluxes, the total thermal fluxes in X and Y directions are 




(17) 


and 
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Introducing the laminar Prandtl number, Pr ,=C p fX / K , into Eqs. 
(17) and (18), the thermal heat flux terms become 


M, Pr, 


Q,=-c, 


Pr, ax n, 


dX J 


where 5 is the Reynolds analogy factor. Introducing Eq. (21) into 
Eqs. (19) and (20), we obtained the following expressions 
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and 


Q,=-c p 


JL?L + ^?L' 

Pr, dX Pr, dX J 


(22) 


The ratio, fl t f fl q1 is usually called the turbulent Prandtl num- 
ber, Pr r In a turbulent pipe flow study [4] it was found that the ratio, 
fi / fl n varied from about 1 at near wall location to about 1.5 at 
the pipe center. Another study [4] of a boundary layer on a heated 
wall also revealed that the ratio, fl q / fi n increased from about 1 at 
the near wall location to approximately 2 at the boundary layer edge. 
Thus, the ratio, fi q / fl n is usually assumed to have a value of 1 at 
near wall locations. The above results were obtained from incom- 
pressible turbulent flow investigations. 

An experiment [5] was recently performed to study the heat trans- 
fer of a shock wave/turbulent boundary layer interaction flow. The 
interaction flow field generated by a 16° compression ramp in a 
Mach 2.84 flow was considered. The results (Fig. 2) show the 
Reynolds analogy factor (the ratio of the Stanton number and half 
of the skin friction factor) increased, significantly at the downstream 
locations. The factor appeared to level-off at a value of 2.15 at the 
downstream location. This may imply that fi q I fi, may be greater 
than I at near wall locations downstream of the compression flow 
field. 

Based on the above observations, we propose that, for the com- 
pression ramp flow the conventional definitions must be modified, 
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S u = S d = 1 
S u = 1 . S d = 2 

S u = 1,S d = 3.65-(X/5 o -2.5)0 38 

Experiment (ref. 5) 
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(23) 


for the thermal fluxes. We further assumed that 5 only varies along 
the X-direction. As shown in Fig. 2, three different potential anal- 
ogy factor variations, 5 as function of X / <5 C , were considered in 
this study; 

1 . The analogy factor is constant, 5=1. 

2. A discontinuity in 5 occurs at the comer X / 8 0 = X C / 8 0 . S is 
1 at upstream locations and 5 changes to a value of 2 at down- 
stream locations. 

3. A discontinuity in 5 occurs at the comer X / 8 =X C / 8 0 . S is 

1 at the upstream locations [X 1 8 Q < X c / . At the down 

stream locations, 5 varies approximately following the ratio of 
the Stanton number and the skin friction factor from the 
experiment [5]. 

In the present study, the Baldwin-Lomax eddy viscosity formu- 
las and the k — E turbulence models were first used to model the 
turbulent momentum fluxes in the Navier-Stokes equations. Then, 
the above Reynolds analogy, with Pr, = 0.93, was used to model the 
thermal flux in the energy equation. Numerical computations were 
performed to solve the transport equations. Since the surface heat 
transfer measurements were only available for the 1 6° compression 
ramp flow, the heat transfer portion of the present computational 
analysis was also limited to the 16° compression flow field. 
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Figure 2. — Reynolds analogy factor, e = 1 6°. 


COMPUTATIONAL METHODS 

The mass averaged unsteady compressible turbulent Navier- 
Stokes equations and the energy equation were solved numerically 
for the steady state properties within the compression ramp flows. 
These properties are the density, p, the velocity components, U and 
F, and the static temperature, T. The equation of state for an ideal 
gas was used to calculate the static pressure, p. The turbulent eddy 
viscosity, Ji r was calculated, with or without the turbulence quanti- 
ties, k and e, based on the selection of the turbulence modeling tech- 
nique. These properties were then used to compute the Reynolds 
shear stress, x XY , skin friction factor, C ^ and the surface heat trans- 

fer ,G h, 


5 


A time-dependent explicit numerical computational method was 
used. Shang et al. [91 developed such a computational code, using 
MacCormack’s finite difference scheme [14]. This code used the 
standard Baldwin-Lomax eddy viscosity concept to model the tur- 
bulence terms. This code was revised by modifying the Baldwin- 
Lomax eddy viscosity concept and by adding the present k — £ 
turbulence models as described in the section TURBULENCE 
MODELS. The revised code with four models available was used 
for the present computations. The conservation equations, the coor- 
dinate transformations between the physical and the computational 
domains, and the finite difference scheme used in the numerical 
computations were described in a previous report [15]. 

Fully developed turbulent flat plate boundary layer flow proper- 
ties were used as the upstream boundary conditions at XI 8 a = 0 . 
These boundary flow properties were first calculated by using the 
computational code. The turbulent flat plate l/7th power law veloc- 
ity profile was assumed for the initial stream-wise mean velocity 
components. The initial normal direction mean velocity components 
were set to zero. The velocity and temperature relationship of a 
compressible laminar boundary layer along a heated flat plate [15] 
was used to define the initial mean temperature distributions. The 
far-field boundary conditions were prescribed by the free -stream 
conditions. The no-slip condition, constant wall temperature, and 
zero normal pressure gradient were used for the wall boundary con- 
ditions. It was further assumed that the flow properties did not change 
along the external mean velocity direction. With this assumption, 
zero property gradients along the X-direction were used as the up- 
stream and downstream boundary conditions. The physical compu- 
tational domain for establishing the flat plate boundary layer condi- 
tions was 7.5 8 0 in the X-direction and 2 5 0 in the F-direction. An 
//-type grid with a 40 by 149 mesh size was used in the computa- 
tions. The grid points were packed near the wall along the F-direc- 
tion. However, the grid points were uniformly distributed along the 
X-direction. For the k and £ computations, the incompressible flat 
plate boundary layer turbulent kinetic energy [ 1 6] was used as the 
initial conditions of k . The initial dissipation rate was calculated by 
using E-pC^k 2 i p, where Jl f was obtained from the Baldwin- 
Lomax turbulence model computations. 

The computed turbulent flat plate boundary layer flow properties 
were used both as the upstream boundary conditions and the initial 
conditions of the Navier-Stokes computations of the compression 
ramp flows. The physical domain was chosen such that the reflected 
shock wave crossed the downstream out-flow boundary and , there- 
fore, the free-stream properties were used as the far-field boundary 
conditions. With the present free-stream conditions, the 8° and 16° 
compression ramps induced slightly different reflected shock wave 
orientations. Consequently, two different physical domains were 
chosen respectively for the 8° and 16° compression ramp flow com- 
putations. A domain of 7 5^ in the X-direction and 2 8 a in the In- 
direction was used for the 8° compression ramp computation, and a 
domain of 7.5 8 Q in the X-direction and 2.5 8 a in the F-direction was 
used for the 16° compression ramp computations. The flow proper- 
ties at the out-flow boundary were also assumed to remain constant 
along the downstream free stream direction. An H - type grid with a 
221 by 149 mesh size was used for the present computations. 

A general concept [15] of the transformation between the physi- 
cal and the computational domains was also used in the present 


computations. The transformation was required to concentrate the 
computational mesh points within the near wall and the comer re- 
gions where large gradients in the flow properties would occur. Two 
different grid configurations were used previously [8] for the 8° 
and 16° compression ramp flow computations. These grid configu- 
rations were optimized, by changing the packing parameters [15], 
to obtain satisfactory computational results of the Reynolds shear 
stresses. These grid configurations were also used here. Using the 
computed flat plate boundary layer flow properties as the initial 
conditions, the numerical computation, with a CFL number = 0.4, 
required approximately 45,000 time iterations to obtain a steady 
state solution. The computation process was time consuming. 

RESULTS AND DISCUSSIONS 

Based on the turbulence modeling techniques, some of the com- 
putational results of the flow properties, such as the turbulent eddy 
viscosity, the Reynolds shear stress, the stream-wise mean flow 
velocity, and the surface heat transfer rate, are described and dis- 
cussed in the following sections. To assess the accuracy of the present 
computations, the computational results are also compared with simi- 
lar results from experiments [5, 10, and 11]. 

Upstream Flat Plate Boundary Layer Flow Properties 

The computed Reynolds shear stress profiles, at X / 5 0 =0.5 lo- 
cation, with the standard and modified Baldwin-Lomax turbulence 
models are shown in Fig. 3(a). A constant Reynolds analogy factor, 
S = 1, was used at this X-location in the computations. The compu- 
tational results are also compared with existing Reynolds shear 
stresses measurements [10]. Both models predicted approximately 
the Reynolds shear stress variation. The standard Baldwin-Lomax 
model calculated larger Reynolds shear stresses at locations away 
from the wall surface. By changing the value of from 0.3 [6] 
to 0.8, the modified model could predict the experimental Reynolds 
shear stresses at the locations away from the wall surface. Thus, a 
different intermittency factor is required to use the Baldwin-Lomax 
model for supersonic boundary layer Reynolds shear stress compu- 
tation. 

The computed Reynolds shear stress profiles at X / 8 0 — 0.5 from 
the computations with low-Reynolds-number and hybrid k — £ tur- 
bulence models were plotted and compared with the measurements 
in Fig. 3(b). The computational results agree very well with the 
measurements. Particularly noteworthy, these turbulence modeling 
techniques could describe the near wall Reynolds shear stress varia- 
tion as it was shown by the experimental results. 

The present computational results of the boundary layer mean 
velocity profile at X / 8 a = 0.5 were studied, using an existing com- 
pressible boundary layer flow analysis [17]. This existing analysis 
used a transformation theory [18] to map the compressible turbu- 
lent boundary layer flow velocity profile into an incompressible 
flow law of the wall velocity correlation. A brief description of the 
mapping technique was presented in the APPENDIX of the paper. 
This mapping technique was used to correlate the F-directional varia- 
tions of the velocity component, U, and density, p, from the present 
computations. The m apping process leads to the prediction of the 
skin friction factor. This skin friction factor can be compared with 
the skin friction factor obtained from the computations. The experi- 
ments [11] measured the fiat plate boundary layer velocity and Mach 
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Reynolds shear stress, (txyA>«. JxlO 3 



Reynolds shear stress, {xxy/Poo Llf, )x1 0 3 
Figure 3. — Upstream Reynolds shear stress profiles, X/6 0 = 0.5. 

number profiles. Assuming an adiabatic wall condition, we calcu- 
lated the boundary layer density variation based on the experimen- 
tal Mach number profile. The mapping technique was then used to 
correlate the velocity component, U, and the density, p. This corre- 
lation was then compared with the correlation based on the present 
Navier-Stokes computational results of U and p. The law of the 
wall, in terms of the velocity scale, U t U x , and the transformed 
coordinate, £, are shown in Figs. 4(a) and (b). Using the Baldwin- 
Lomax models, the computations predicted larger U I U r than that 
was prescribed by the law of the wall (Fig. 4(a)). Using thej^ - £ 
turbulence model, the computations predicted a smaller U I U T than 
that was prescribed by the law of the wall (Fig. 4(b)). Using the 
hybrid k — E turbulence model, the present computations predicted 
large U / U z at locations near the viscous sublayer edge (Fig. 4(b)) 
and failed to resolve the velocity profile within the buffer layer. 
This may be because the model did not consider the turbulence mod- 
eling within the buffer layer. Despite this, the hybrid k - £ model 
generally does a good job of describing the flow properties. A value, 
Cy =0.002 , for the incompressible skin friction factor was required 
to establish the correlations from the present computational results 
of the velocity component, (/, and the density, p. Based on this C y 
value, Eq. (A5) gives C f =0.001 . The present computations also 
predicted a skin friction factor, C f =0.0009. In addition, the com- 



Transformed coordinate, £ 



Figure 4. — Law of the wall velocity profiles, x/5 0 = 0.5. 

putations calculated a surface heat transfer rate, 

G =0.75 Btu/ft 2 - sec, for a near adiabatic wall condition, 

w 

T w / T 0 „= 1.04. The above observations indicated that the present 
Navier-Stokes computations could predict the surface heat transfer 
and the skin friction factor. 

As it was shown by the results in Figs. 4(a) and (b), the experi- 
mental data agrees well with the law of the wall for £<10.6 • F° r 
£> 1 0.6*. the experiment shows a larger velocity scale than that 
was calculated by using Eq. (Al). With the Baldwin-Lomax model 
for turbulence modeling, the computational velocity (Fig. 4(a)) 
agrees closely with the experimental data. With the k — E turbu- 
lence transport equations for turbulence modeling, the computa- 
tional velocity correlations (Fig. 4(b)) show a lower velocity scale 
than that was indicated by the experimental data. At large £ loca- 
tions, the computations also predicted a lower U / U T than that was 
given by Eq. (Al). 

The present computational results of the eddy viscosity from dif- 
ferent turbulence models were also plotted in terms of the trans- 
formed coordinate £ in Figs. 5(a) and (b). The modified Baldwin- 
Lomax model (with C klcb “ 0 .8) calculates small eddy viscosity at 
the location away from the surface (Fig. 5(a)) and this model pre- 
dicts the experimental measurements of the Reynolds shear stress 
(Fig. 3(a)). The low-Reynolds-number k - £ turbulence model and 
the hybrid k — E model calculate approximately the same eddy vis- 
cosities within most parts of the boundary layer (Fig. 5(b)). A dif- 
ference in the eddy viscosities occurs only at locations very close to 
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Eddy viscosity, iVp Ec *dy viscosity, iVp 


Baldwin-Lomax model, S = 1 



Transformed coordinate, £ 



Transformed coordinate, £ 

Figure 5. — Upstream eddy viscosity profiles, x/5 0 = 0.5. 


the surface. These k — € models only induce differences in the 
Reynolds shear stress predictions at the near surface locations (Fig, 
3(b)). A comparison between the results of eddy viscosity in Figs. 
5(a) and (b) also reveals that both the k - E models predicted ap- 
proximately die eddy viscosity as it was calculated by using the 
modified Baldwin-Lomax eddy viscosity model. This also verifies 
that the Baldwin-Lomax eddy viscosity formulation requires 
Cjjeb = 0*8 in order to model the turbulent momentum flux of a 
supersonic boundary layer flow. 

8 Degree Compression Ramp Downstream 
Flow Properties 

The computed Reynolds shear stress profiles, based on the stan- 
dard and modified Baldwin-Lomax turbulence models and the hy- 
brid k-t models, at five different X 1 8 o locations were shown in 
Fig. 6. A constant Reynolds analogy factor, S u =S d = 1, was used 
in the computations. The computational results were also compared 
with the experimental measurements [10]. The standard Baldwin- 
Lomax model predicted lower downstream x XY / p^U* values than 
were measured. With C i = 0.45 and C kJeb - 0.8, the modified model 
could predict accurately the Reynolds shear stress distributions 
downstream of the shock wave. The present computations also re- 


quired a larger C- value to account for the downstream Reynolds 
shear stress amplification. Based on the hybrid k — E turbulence 
modeling techniques, the computed Reynolds shear stresses also 
agreed well with their measurements especially at the downstream 
locations. Both the computation and the experiment gave down- 
stream Reynolds shear stress profiles which are significantly dif- 
ferent from the profile at upstream location. Along the F-direction 
and near the shock wave locations, the computed Reynolds shear 
stress increases and then decreases to its free stream condition. 
However, the measurements indicate that, near the shock wave lo- 
cations, reduces and then increases to its free stream value. 

The computed surface pressure, p w f p„, and the skin friction 
factor, Cj = T W / (1 / 2 )p ao Ul, variations are shown in Fig. 7. The 


surface shear stress was computed by using 



velocity components U and V were used to obtain numerically the 


normal velocity gradient, ^ 


, at the wall surface. At an upstream 


location (X < X c ), the velocity gradient, 



dU 

dY 


was evalu- 


ated based on the U values at the first three grid points from the 
surface. At a downstream location (X > X c \ V and V at the first 
three grid points were first used to determine the velocity compo- 
nents U s at the grid points. These velocity components were then 

used to evaluate numerically the normal velocity gradient, ^ 


The experimental results [10] of these surface conditions were also 
plotted in Fig. 1 0 for comparison. Using either Baldwin-Lomax tur- 
bulence model, the computations predicted the experimental sur- 
face pressure distribution. Using the modified Baldwin-Lomax 
model, the computations calculated a higher level of the skin fric- 
tion downstream of the compression comer. This skin friction dis- 
tribution agrees better with the measurements. Using the hybrid 
k — z model, the computations could also predict accurately both 
the experimental surface pressure and skin friction variations. 


16 Degree Compression Ramp Downstream 
Flow Properties 

Using the best turbulence models for the 8° ramp and the Reynolds 
analogy assumptions, the computed Reynolds shear stress varia- 
tions along the F-direction at five different X i 5 o locations were 
shown in Figs. 8(a) and (b). For comparisons, the corresponding 
shear stress measurements [10] were also plotted in Figs. 8(a) and 
(b). Figure 8(a) shows these results when the modified Baldwin- 
Lomax model and different Reynolds analogy assumptions were 
used respectively in the computations. With constant Reynolds anal- 
ogy factor, S u = S d = 1, the computations well predicted the experi- 
mental Reynolds shear stress variations along the F-direction at 
X f 3 0 = 5.5 and 6.5 locations. However, the computed Reynolds 
shear stresses, at the near comer locations, were larger then the 
measured Reynolds shear stresses. With the experimental Reynolds 
analogy factor, S d =3.65- (X / 8 o -2.5) , the computations 

overpredicted the Reynolds shear stresses at almost all downstream 
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Baldwin- Lomax model, S = 1 

Modified Baldwin-Lomax model, S = 1 

Hybrid k-e model, S = 1 

O Experiment (ref. 10) 



Figure 6.— Reynolds shear stress profiles, 8 = 8°. 



Figure 7. — Surface pressure and skin friction factor of 8 = 8°. 


locations. Figure 8(b) presented the results when the hybrid k - £ 
turbulence model and different Reynolds analogy assumptions were 
used in the computations. The computations, using S u = S d - 1, 
computed approximately the experimental T-directional shear stress 
profile at X / 5 0 =3 location. Using either S d = 1 or S d = 2 in the 
computations, the computations could predict the shear stress varia- 
tion along the T-direction at X 1 5 o =6.5 location. Changing the 5^ 
value from 1 to 2 in the computations, the computations produced 
small differences in the Reynolds shear stress predictions. 


The computed skin friction factor and surface pressure variations, 
with different turbulence models and the Reynolds analogy factor 
assumptions, were plotted and compared with their experimental 
measurements [10] in Figs. 9(a) and (b). The computational results 
of the skin friction factor were obtained in a manner as it was used 
in the 8° compression ramp flow calculations. These results indi- 
cated that the downstream Reynolds analogy assumption did not 
significantly affect the present computations of the surface pres- 
sure and skin friction factor. Using the modified Baldwin-Lomax 
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Modified Baldwin-Lomax model, S u = = 1 


— Modified Baldwin-Lomax model, S u = 1 , = 3.65-PC/80-2.5) 0 - 38 

* Experiment (ref. 10) 



— Hybrid k-e model, S u = 1 , = 2 

Hybrid k-e model, S u = 1 , = 1 

O Experiment (ref. 10) 



model (Fig. 9(a)), the computations closely predicted the surface 
pressure but poorly predicted the skin friction factors. The compu- 
tational results of the skin friction factor indicated that flow separa- 
tion occurred near the comer location. Just the opposite was true, 
using the hybrid k — Z turbulence model (Fig. 9(b)), the computa- 
tions calculated a lower level of the downstream surface pressure 
and accurately predicted the skin friction factor variation. The com- 
putational results of the skin friction factor also indicated that a 
very small flow separation region existed around the comer loca- 
tion. This is close to the experimental observation [101 that incipi- 
ent flow separation could occur near the comer location. 

For a compression ramp flow, it has been shown [19] that the 
upstream and downstream static pressure ratio played an important 
role in modifying the turbulence within the flow field. The role of 
the pressure ratio was not included in the hybrid k — E turbulence 


model. It could not capture fully the turbulence distortion in the 
compression ramp flow and the present computation could not pre- 
dicted accurately the surface pressure on the 16°compression ramp. 
The success of the present calculation of the skin friction along the 
16° ramp may require further investigation. However, the hybrid 
k — Z turbulence model is straightforward and is a worthwhile ad- 
vancement for the engineering application. 

The Navier-Stokes computational results of the static tempera- 
tures at the first three grid points from the wall surface were used to 


at the 


dr 

calculate numerically the normal temperature gradient, 
wall surface. This temperature gradient was then used to compute 

n --K ?L 

the surface heat transfer rate by using fa ' The com ” 


f 
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Experiment 
(ref. 10) 

□ Surface pressure 

O Skin friction factor 

Modified Baldwin- Lomax model 


S u = S d = 1 

S u = 1 l S d = 3.65-(X/5 0 -2‘5)°- 38 



Stream-wise distance, X/5 0 

Experiment 
(ref. 10) 

□ Surface pressure 

O Skin friction factor 


Hybrid k-c model 
S u = 1 , S d = 1 



St ream -wise distance, X/8 0 

Figure 9. — Surface pressure and skin friction factor, 8 = 16°. 


puted surface heat transfer results, based on different assumptions 
of the Reynolds analogy factor and turbulence models, were plot- 
ted and compared with the measurements [5] in Fig. 10. With con- 
stant Reynolds analogy factor, 5 = 1, the computational results did 
not show significant increase in the surface heat transfer at the down- 
stream locations. With the experimental Reynolds analogy factor 
and the modified Baldwin-Lomax turbulence model, the computed 
surface heat transfer rate decreased at locations just upstream of the 
comer location. The surface heat transfer rate oscillated at loca- 
tions immediately downstream of the comer and, then, increased 
rapidly to a higher value than that was measured. Using the hybrid 
fc — e turbulence model and S u = S d = 1 assumption, the computa- 
tions even predicted a reduction in the surface heat transfer at the 
downstream locations. However, using a S d = 2 assumption, the 
computation predicted the downstream surface heat transfer varia- 
tion. A turbulent Prandtl number Pr t = 0.93 was used in the compu- 


Hybrid k-c model 

S u = 1 , S d = 1 

S u = 1,S d = 2 


Modified Baldwin-Lomax model 




Figure 1 1 . — Correlations of velocity and temperature fluctuations, 
8 = 16 °. 


tations. Thus, the combination of the hybrid k — £ turbulence model 
and a ratio of JJ q f fi t =2.\5 seems to be an appropriate turbulence 
modeling technique for the computations of the surface heat trans- 
fer rate downstream of the 16° compression ramp. 

It should be noted that for one case, the case with S u = 1 and 
S d =3.65- (X / 5 a -2.5) assumptions and the modified 
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Baldwin-Lomax turbulence model for 16° ramp, only preliminary 
computational results are described in this paper. Because of the 
very long computing time, we used the previously computed 8° 
compression ramp mean flow properties to start this part of compu- 
tations. The computational results did not show significant varia- 
tions after 15,000 time iterations. Thus, the computational results, 
at the 15,000th iteration, were presented in Figs. 8(a), 9(a), and 10. 

Equations (15), (16), and (21) were used to model the turbulent 
heat flux terms for the present computational analysis. The 16° com- 
pression ramp computations (with S u = 1 and S d = 2) indicated the 
turbulent heat flux /p f JJT changed significantly as flow moved 
downstream. For example, the computed profiles of the v 't ' / UT 
correlation at the upstream (X/ 8 o = 0.5) and the downstream 
X l 8 a =6.5 locations were plotted in Fig. 11. These profiles indi- 
cated that, within the boundary layer flow, the correlation v 'f / UT 
was amplified at the downstream location. 

Recently, advanced turbulence modeling techniques have been 
developed to model the near wall turbulence. Turbulence transport 
equations [20] were used to analyze the near wall turbulent heat 
flux. The turbulence modeling techniques could calculate the ve- 
locity-temperature correlations within some incompressible bound- 
ary layer flows. An analysis [21] was also proposed to relate the 
evolution of the Reynolds stress and the turbulent heat flux to the 
mean flow properties in a shock- wave/turbulence interaction. These 
new turbulence modeling techniques may lead to a better analytical 
approach for the predictions of the Reynolds shear stress, the turbu- 
lent heat flux, and the surface heat transfer rate of a shock-wave/ 
turbulent boundary layer interaction flow field. 


(c) The ratio of the transport coefficients used in the Reynolds 
analogy has an impact on the computation of the surface heat trans- 
fer rate downstream of the compression ramp. With a transport co- 
efficient ratio of 2/Pr, and a hybrid k - £ turbulent momentum 
flux modeling technique, the present computations predicted the 
experimental results of the heat transfer rate on a 16° compression 
ramp surface. 


APPENDIX 

Baronti and Libby [17] performed an analysis to map a com- 
pressible turbulent boundary layer flow into an incompressible 
boundary layer flow. Their results showed that the compressible 
flow velocity profile was well correlated by the law of the wall. The 
law of the wall was written as 


U/U x =& for 0<C<C/ 
= 2.43 In 7.5 for l f <£ 


(Al) 


where £/ t =(t„ /p w )° 5 is frictional velocity. f / =10.6 was as- 
sumed. The y-directional transformed coordinate, £, was defined 
as 



CONCLUSIONS 

The application of turbulence modeling techniques and the 
Reynolds analogy in the Navier-Stokes computations of two-dimen- 
sional supersonic compression ramp flows were investigated. The 
Baldwin-Lomax eddy viscosity concept and the k — £ turbulence 
models (with and without modifications) were used for turbulent 
momentum flux modeling. The Reynolds analogy was used for the 
turbulent thermal flux modeling. The ratio of the thermal and mo- 
mentum transport coefficients was assumed to be Sf Pr r S is an anal- 
ogy factor and it was also assumed to be a function of the stream- 
wise location. The turbulent Prandtl number is constant, Pr, = 0.93. 
The Navier-Stokes equations and the energy equation were solved 
numerically for the flow properties of the compression ramp flow 
field. The Reynolds shear stress, skin friction factor, and the sur- 
face heat transfer rate were calculated and compared with their mea- 
surements. It was concluded that: 

(a) The Navier-Stokes computations predicted the supersonic flat 
plate boundary layer flow properties when the k and £ transport 
equations were used for turbulence modeling of a boundary layer 
flow with nonadiabatic wall condition. 

(b) A modified Baldwin-Lomax turbulence model and a hybrid 
£ _ £ turbulence model were found to be effective for the Navier- 
Stokes computations of the Reynolds shear stress and the skin fric- 
tion of a supersonic compression ramp flow. The ratio of the trans- 
port coefficients in the Reynolds analogy only slightly influenced 
the computations of the Reynolds shear stress, the skin friction fac- 
tor, and the surface pressure. 


with 


M 

cm 

ii 

b 

s 

M 

of 

f p ~] 

^ P ) KP-) 

U/J 

J Jo ! 

P J 




(A3) 


Some relationships between the compressible and incompress- 
ible flow properties were also written as 


and 



(A4) 




p./O 

/ \ 

V 

(p.pj 

Uw-J 




(A5) 


Given a Cj value and the compressible flow velocity and density 
profiles, Eqs. (A2) to (A4) could be used to calculate U / U T and_£ 
at a specific Y location. It is also possible to select a value of C f 
which correlates U / U x and £ in the law of the wall region. Based 
on the C f value, Eq. (A5) could then be used to calculated com- 
pressible flow skin friction factor 
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